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We show how recent angle resolved photoemission measurements addressing the Fermi arcs in the cuprates 
reveal a very natural phenomenological description of the complex superfiuid phase. Importantly, this phe- 
nomenology is consistent with a previously presented microscopic theory. By distinguishing the order parame- 
ter and the excitation gap, we are able to demonstrate how the collapse of the arcs below T c into well defined 
nodes is associated with the smooth emergence of superconducting coherence. Comparison of this theory with 
experiment shows good semi-quantitative agreement. 

PACS numbers: PACS numbers: 74.20.-z, 74.20.Fg, 74.25.Bt, 74.25.Fy 
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In the past most of the interest in lower T c cuprate super- 
conductors has focused on the exotic, non-Fermi liquid nor- 
mal phase. In a recent paper by Kanigel et al, 1 angle resolved 
photoemission spectroscopy (ARPES) experiments on several 
underdoped samples of I^S^CaC^Os (Bi2212) were used 
to establish key features of the superconducting phase. In par- 
ticular, it was reported 1 that (i) the ARPES -measured exci- 
tation gap, A(k,T) is roughly constant in temperature from 
T = to above T c . (ii) Below T c , A(k) displays the d-wave 
point nodes which broaden into Fermi arcs above T c , with the 
change occurring within the width of the resistive transition at 
T c . (iii) It is claimed 2 that the energy scale of the excitation 
gap is T*, or the pseudogap 3,4 onset temperature, and that the 
Fermi arc length scales with T/T* above T c . From (i) it is 
inferred that (iv) "the energy gap is not directly related to the 
superconducting order parameter". 

These latest experiments have underlined the fact that the 
superfiuid phase is itself very complex in the presence of a 
normal state gap or "pseudogap". 3,4 Some theories 5 - 6 seem to 
suggest that the already large pseudogap becomes the order 
parameter immediately below T c , which would seem to im- 
ply an (unphysical) jump in the order parameter and in the 
superfiuid density, n s . In our approach we show how these 
important photoemission observations reveal a more natural 
description of the superfiuid phase. We then review a micro- 
scopic model consistent with this phenomenology which has 
been demonstrated 4 to be compatible with a variety of other 
experiments and show that it yields good semi-quantitative 
agreement with a large number of different representations of 
these recent photoemission data. 

Our microscopic scheme is based on a BCS-Bose-Einstein 
condensation (BEC) crossover scenario 4,7 and is distinct from 
the phase fluctuation scenario. 5,8,9 It has been argued 10 to be 
appropriate to the cuprates because of their very short co- 
herence length. We emphasize the generic features of our 
framework. One assumes that there are attractive interactions 
which lead to pairing which, in turn, gives rise to a gap in the 
fermionic spectrum. Noncondensed pairs above T c contribute 
to this excitation gap, A, just as do condensed pairs. Thus, we 
infer that below T c , A contains two contributions, A sc from 
condensed and A pg from noncondensed pairs. Since the pair 
density is associated with the square of the gap, these contri- 



butions add in quadrature to yield 
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The same equation is consistent with points (i) and (iv) above. 
Note that the superfiuid density, n s (T) is observed to van- 
ish smoothly as T c is approached from below, which requires 
that the superconducting order parameter A sc turn on continu- 
ously as in a second order phase transition. We then conclude 
that because the excitation gap A(k, T) is roughly a constant 
in T across T c there must be another component to the ex- 
citation gap below T c , which compensates for the T depen- 
dence in A SC (T). The simplest approach is to think of this 
term (i.e., A 2 ff (T)) as a "fluctuation" contribution of the form 
(A 2 ) - (A) 2 . 

The normal state analysis of ARPES experiments has al- 
ready made substantial us o 6 ' 11 ' 12 of a broadened BCS form 13 
for the fermion self energy 
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Here the broadening 7 ^ and "incoherent" background 
contribution E refect the fact that noncondensed pairs do 
not lead to true off-diagonal long-range order. We adopt a 
tight binding model for the band dispersion e^, although the 
detailed band structure, is of no importance in our calcula- 
tions which address ARPES data points along the Fermi sur- 
face ek = 0. We define Ak lPg = A pff (^k, and we introduce 
</3k = cos(20), to reflect the d-wave k dependence along the 
Fermi surface. 

In contrast to earlier scenarios 5,6 here we distinguish the su- 
perfiuid from the normal phases via an additional component 
of the self energy arising from the condensate, E sc so that 



E(k,w) = E pg (k, lu) + E sc (k, uj) 



E sc (k,w) = 



k.sc 



(3) 
(4) 



with Ak. sc. = A sc <,5k- E sc is associated with long-lived con- 
densed Cooper pairs and so it is of the same form as E pg 
but without the broadening. With this self energy in the 
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Green's functions, the resulting spectral function, A(k, ui) = 
— 2ImG(k, u> + iO) can be readily determined. One can see 
that the spectral function at all k contains a zero at u> = — £k 
below T c , whereas it has no zero above T c . This dramatic 
effect of superconducting coherence will be reflected in the 
disappearance of Fermi arcs below T c . 

Since A pg represents a contribution from thermal bosonic 
fluctuations, it should vanish in the ground state and roughly 
be of the general form 



K q {T) = (T/T c )"A 2 (T c ), T<T C 



(5) 



Using this equation one can determine T c as the the lowest 
temperature at which A sc = 0. The exponent v in general 
varies between 3 and 3/2, which corresponds to linear and 
quadratic pair dispersion, respectively. What is important is 
not the value of v, but that at T c the order parameter vanishes, 
and the term A pg accounts for the entire excitation gap A(T C ). 

In the underdoped regime, because A(T) « A(T C ), it 
would appear that the smooth vanishing of the superfluid den- 
sity n a at T c introduces a challenge for theory. In one-gap 
scenarios 5,6 [which seem to assume A SC (T~) = A(T~)] one 
needs to demonstrate explicitly why there is not a disconti- 
nuity in n s at T c . By including two gap parameters below 
T c , we arrive at a natural understanding of n s . Importantly, 
it is only the condensed pairs which contribute to the Meiss- 
ner magnetic screening, as expected, so that n s oc A 2 C , which 
thus vanishes smoothly at T c . Indeed, this is a natural conse- 
quence if one considers superconductivity as Bose condensa- 
tion of Cooper pairs. The noncondensed pairs discussed above 
contribute to the destruction of n s in addition to the usual 
fermionic terms. 14 These bosonic contributions are clearly re- 
quired since the number of fermionic excitations (oc T/A) is 
smaller than the total fermion number for all T < T c . The 
phase fluctuation scenario 8 similarly invokes bosonic excita- 
tions to destroy n s at T c . 

We will show below that, at a semi-quantitative level, we 
can address ARPES data by assuming for all T <T* 
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Here E k = V^k+^k' with = A^ k ; U < is the 
pairing interaction strength, 14 and f(x) is the Fermi function. 
The ratio A(T) /A(0) as a function of T/T* is then given by 
a nearly universal curve independent of the doping concen- 
tration, x. In the underdoped regime, A(T C ) is large so that 
Eq. (|6]l is consistent with a very high T* » T c and a roughly 
constant gap within the superfluid phase. 

We briefly present a T-matrix approach from which 
Eqs. dU - © were previously derived, noting that more 
details can be found elsewhere. 4 The BCS-like constraint in 
Eq. d6) can be interpreted as equivalent to the BEC condi- 
tion that the noncondensed pairs have zero chemical potential 
Impair = 0. at and below T c . This determines the form of the 
noncondensed pair propagator t pg (Q) = U/[l + Ux(Q)] such 
that t P9 (0) = oo at T < T c , where X (Q) = J2k G o(Q ~ 
K)G(K)(p 2 ,_ ci ^ 2 , is the pair susceptibility. Here G and Go 
are the full and bare Green's functions, respectively. To make 
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Figure 1: (color online) Spectral function A(u>) for x = 0.125 at 
(a) T — T c and (b) 0.95T C at different angles cj) along the Fermi 
surface of a cuprate superconductor, where <f> increases by 4.5° from 
to 45° along the direction of the arrow. The upper inset shows the 
temperature dependence of A(uj) at <f> = 32° near the node, with 
T/T c = 1, 0.99, 0.95, 0.9, 0.8, 0.6 and 0.3 in the direction of the 
arrow. Shown in the lower inset are the gaps vs T below T c . The 
effects of phase coherence are more pronounced in the nodal region 
of the Fermi surface, where the gap is small. 



direct association with Eq. @, we drop the 7 and So term 
(which are more important above than below T c ) in Eq. (O 
and arrive at E(k, J) « A£/(w + e k ) for T < T c . This rea- 
sonable approximation can be used to determine the form for 
G, and establish an equivalence between Eq. (|6]l and the BEC 
condition fi pa ir = 0. 

To establish the validity of Eqs. (O and dU, we note that 



there are two contributions to the full T-matrix t = t 



pa 



t. 



A 2 

pg 

A 2 



where the condensate contribution t sc (Q) — — 7fr L 5(Q). 
Similarly, we have for the fermion self energy 'E(K) = 
V SC {K) + Tipg(K) = J2 Q t(Q)Go(Q - K)yl_^, from 
which Eq. (01 follows at once. A vanishing chemical po- 
tential means that t pg (Q) diverges at Q = when T < 
T c . Thus, in the superfluid phase only, we approximat e) 13 ' 15 
£„„(#) to yield ^ pg {K) » -G (-K)Al pg , where 
= - Y^Q^pgiQ)- We thus may write E ps (k,o;) = 
k pg(T)/(w + £k) + small corrections, where we accommo- 
date the corrections 13 with the broadening factor 7 and ad- 
ditional term So- This is necessary in order to address the 
concrete fermion spectral function. In this way, Eq. (0 then 
follows. Finally, at small four-vector Q = (ft, q), we expand 
the inverse of t pg , after analytical continuation, to obtain a 
simple quadratic dispersion 4 for the pairs Sl q ks q 2 /(2M*), 
implying v = 3/2. 

To illustrate the simple physics, we do not attempt to do 
detailed curve fitting. We use one set of parameters for all 
doping x. Thus, in our numerical results we take only Eq 
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and -f(95K) as adjustable to optimize overall fits to the mul- 
tiple data sets. More concretely, we obtain the universal curve 
A(T) / A(0) as a function of T/T* from Eq. ©, say, at opti- 
mal doping, and use the experimentally known value of A(0) 
at given doping concentration x to determine an input T* 
and input gap A(T). This value of T* is to be distinguished 
from the T* obtained in ARPES experiments, here called T e * . 
For our tight binding, quasi two-dimensional lattice, Eq. (_6]i 
yields 2A(0)/T* w 4.3. In ARPES experiments, T* x is de- 
termined as the temperature where 7 and A are roughly equal 
so that there is no observable density depletion at the Fermi 
level. The experimentally deduced ratio from ARPES data 2 
is slightly larger than 5 for moderately underdoped samples, 
and increases with underdoping. A higher ratio of 8 has been 
observed in local STM measurements. 16 The next step the- 
oretically is to use the known T c to determine A pg (T) and 
A SC (T) below T c , using Eqs. (0 and (___}. One can then de- 
termine an "experimental" T e * and spectral gap following the 
ARPES prescription for given parameters 7 and S . 

To make this determination, we choose our adjustable pa- 
rameters £0 and 7 via a rough fit to ARPES data near T c (« 
95K) at optimal doping. Since Eq is primarily governed by 
the particle-hole channel, 11 we take it to be independent of 
doping and T, and given by S = 26 meV. The broaden- 
ing parameter 7 depends primarily on temperature. Consis- 
tent with scattering rate measurements in the literaturer^Ui 
we take 7 = 26 meV at 95 K, with 7 = 7 (95K)(T/95K) 
above T c and 7 = ^(T C )(T/T C ) 3 below T c for given doping 
x, as used earlier. 19 . Finally, this yields a ratio 2A(0)/T e * = 
5.6 ~ 5.7 for the cases presented below, consistent with 
experiment. 2 Throughout = and 7r/4 = 45° denote the 
anti-nodal and nodal directions, respectively. In order to com- 
pare with ARPES data, we convolve the spectral function with 
a Gaussian instrumental broadening curve with a standard de- 
viation cr = 3 meV, given by the ideal ARPES resolution. 

In the inset to Fig. 02b) we plot the various gap parameters 
A sc , A and A pg as a function of temperature below T c and 
slightly above, for a typical underdoped system. Here for def- 
initeness we have chosen v = 3/2 in Eq. (0. It can be seen 
that, as T is lowered below T c , the total excitation gap A is es- 
sentially a constant while A sc increases from at T c to reach 
the full gap value at T = 0. By contrast, A pg monotonically 
decreases to 0. 

In Figs. Q~f a) and[TJb), we plot the spectral function A(ui) 
for £k = (on the Fermi surface) at and slightly below T c , re- 
spectively, and for different angles. To illustrate the physics, 
here we do not include the ARPES instrumental broadening. 
Just below T c , the sharp dip at uj = is associated with the 
onset of a very small condensate, which nevertheless leads 
to a depletion of the spectral weight at the Fermi level. The 
coherence associated with the order parameter is better illus- 
trated near the nodal region where a gap, absent at T c , appears 
as T decreases, as shown for <f> — 32° in the upper inset. A 
closer inspection of the shape of A(lo) in Fig. |TJb) suggests 
that A(u>) cannot be described by a simple broadened BCS 
spectral function. 

Figure _2]presents a plot of the spectroscopic gap for an un- 
derdoped cuprate with x = 0.125 as a function of angle and 
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Figure 2: (color online) Spectral gap as measured in ARPES as a 
function of angle cj> along the Fermi surface of an underdoped cuprate 
superconductor of doping x = 0.125, at different T/T c between 
0.5 and 1.5. The spectral function has been broadened with a small 
ARPES instrumental resolution of 3 meV. A Fermi arc appears as 
the extended range of zero gap value. The fast departure of the T = 
0.99T C curve from the T = T c curve indicates a collapse of the 
Fermi arc. 
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Figure 3: (color online) Loss of zero energy spectral weight as a 
function of temperature at the antinode, for different hole concentra- 
tions. Upper right inset shows rescaled T dependence, while lower 
left inset indicates <j> dependence for various T. 



for various temperatures below (dashed curves) and above T c 
(solid curves). Here the spectral gap is given by half the peak- 
to-peak distance in the spectral function, and is smaller than 
the input gap A(T). For T < T c , since the gap is roughly T 
independent, the various curves tend to coalesce. Above T c , 
the extended range of zero gap value around = 45° gives 
rise to the Fermi arcs. 612 This is due to the presence of 7. The 
rapid deviation of the (red dashed) T = 0.99T C curve from 
the (black solid) T = T c curve indicates a collapse of the 
Fermi arc, leading to the protected nodes below T c , reflecting 
the emergence of superconducting order below T c . This figure 
should be compared with Fig. 2 of Ref.Q]. 
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Figure 4: (color online) Fermi arc length as a function of T /T c * for 
doping concentrations from optimal to underdoping for a cuprate su- 
perconductor. Fermi arc length is typically finite above T c and drops 
to zero upon the onset of phase coherence. The normal state portions 
of the curves is close to universal, in agreement with Ref. y. The 
comparison in the inset between the theory with a slightly (15%) 
enlarged T c * and experimental data (symbols)^ shows a good semi- 
quantitative agreement. 

In Fig.[3]we plot the relative loss of spectral weight L(<f>) — 
1 — I(0)/I(A<p), defined following Ref. 1, at the anti-node 
<f> = as a function of T for different doping concentrations, 
where I(u>) = A(oj). This determines T* x as where L(0) van- 
ishes. The lower left inset presents plots of L(<fi) as a function 
of 4> for various temperatures while the upper right inset illus- 
trate the anti-node behavior L{4> — 0) as a function of T/T* x . 
Above T c the dependence is linear in T reflecting a linear T 
dependence in 7. Illustrated in the upper inset is a consid- 
erable universality above T c as a function of T/T* x , which 



provides a prediction for future ARPES data analysis. As the 
temperature decreases, the abrupt jump at T c reflects the on- 
set of phase coherence, as in Fig. |2] Here the x = 0.1 and 
x = 0.135 lines are close to the data for the 67K and 80K 
samples in the upper panel of Fig. 4 in Ref.Qj. 

Finally, Fig. [4] shows the sharp collapse of the Fermi arcs 
from above to below T c ; we plot the percentage of arc length 
as a function of T/T e * and for different doping concentrations 
from the optimal to the underdoped regime. The small non- 
vanishing arc length at low T reflects the finite ARPES res- 
olution. As a consequence of the Fermi arc collapse below 
T c , the nodes are "protected". In addition, there is a clear uni- 
versality seen in the normal state, in good agreement with the 
central features observed experimentally, shown in Fig. 4 of 
Refs. [l| and |2j Since the mean-field equation is only a crude 
approximation above T c , it is reasonable to allow T e * to vary 
slightly, as we have done here, to compare semi-quantitatively 
with the data. By contrast with Ref. 1, however, our curves are 
not strictly straight lines, reflecting the nonlinearity of the gap 
as a function of cf> — 7r/4, as also found in Ref. 6; this seems 
to be consistent with Fig. 4 of Ref.0. 

The microscopic approach 4 considered here is associated 
with stronger-than-BCS attractive interactions which lead to 
small pair size. A major consequence of this theory (as well 
as the equivalent phenomenology which ARPES experiments 
lead us to infer) is that pseudogap effects persist below T c in 
the form of noncondensed pair excitations of the condensate. 
Importantly, this leads to two contributions to the self energy 
below T c [see Eq. ©]. We find that the collapse of the Fermi 
arcs is not to be associated with an abrupt disappearance (as 
assumed elsewhere 6 ) of the inverse pair lifetime 7, appearing 
in S pg , but rather it reflects the gradual emergence of the con- 
densate, appearing in E sc , to which the finite momentum pairs 
are continuously converted as T decreases. 

This work was supported by NSF Grants No. PHY- 
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